NASA TM X-3108 


NASA TECHNICAL 

MEMORANDUM 





CONTROL SYSTEM DESIGN 

USING FREQUENCY DOMAIN MODELS 

and parameter optimization, 

WITH APPLICATION TO 
SUPERSONIC INLET CONTROLS 


by Robert C. Seidel <tnd Bruce Lebtinen 

Lewis Research Center 
Cleveland, Ohio 44135 



HATIONAl AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, 0. C. . OCTOBER 1974 




l H. |. hi r,., 

N A_SA_TM_X-j)H)H 

2 i mvi'f imii'Mi Arrr’v.mii N«i 

3 i: rijiitmtS Cittifn.'i No 

« I'll'- .Ihil *tl. 

t’ONTROI, SYSTK.M 1)1*: SIGN USING FRMQUKNCY DOMAIN 
MOUNDS AND PAUAMFTKR OPT1. "-'ATION. WITH 

. — p p U i‘ A;noN_;ro supfksonic ini, ft contijods 

6 H,|w)n ;w 

ooTc.iii-.T vr;> i 

0 Ptiforihinij Qfffim/.ition Code 

) Author (s) 

Robert C. Seidel and Bruce Lehllnen 

8 Performing Organt/mioii H*:port No. 

K -7710 

9 Per forming Oitj;ini/ation Name and Addrosj 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio t4135 

10 Wo* Unit No. 

501-24 

11, Contact or Graft! No 

12, Sponsoring Agency Nomr and Address 

National Aeronautics and Space Administration 
Washington. D. • . 20546 

13, Type of Report and Period Covered 

Techuieul Memorandum 

1 A. Sponsoring Agency Code 

15. Supplementary NcvC» 



1G Abstract 

A technique is described for designing feedback control systems using frequency domain 
models, a quadratic cost function, and a parameter optimization computer program. 
FORTRAN listings for tim computer program are included in the report. The approach 
is applied to the design of st.ock position controllers for a supersonic inlet. Considered 
arc a deterministic or random system disturbance and the presence of random measure- 
ment noise. The cost function minimized is formulated in the time domain, but the prob- 
lem solution is Obtained us'ng a lrequcncy domain system description. A scaled and con- 
strained conjugate gradient algorithm Is used for the minimization, In applying the approach 
to a typical supersonic inlet, both optimal proportional -plus -integral (PI) and proportional- 
plus -integral -plus-derivative (PID) i onlrollcrs were calculated. For the inlet considered, 
a single-loop PI controller was iudged to be the most desirable of the various designs con 
sidered. 
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CONTROL SYSTEM DESIGN USING FREQUENCY DOMAIN MODELS 
AND PARAMETER OPTIMIZATION, WITH APPLICATION 
TO SUPERSONIC INLET CONTROLS 
by Robert C. Seidel and Bruce Lehtinen 
Lewis Research Center 

SUMMARY 

This report described a technique for designing feedback control systems using fre- 
quency domain models, a quadratic performance index, and a parameter optimization 
computer program. The approach is applied to the design of a terminal shock position 
controller for a mixed-compression supersonic inlet. The computer program described 
can be used to design controllers for any two-loop linear feedback system having a single 
control input and a cascade structure. The problem formulation assumes the system is 
acted upon by a single deterministic or random disturbance plus random measurement 
noise introduced in each loop. The quadratic performance index chosen to be minimized 
is a weighted sum of averaged square system output, output rate, control, and control 
rate. The performance index is expressed in the frequency domain and is minimized, 
given a frequency domain system description, using a scaled and constrained conjugate 
gradient search algorithm. 

For the supersonic inlet problem, the disturbance is a deterministic airflow pertur- 
bation at the diffuser exit, and measurement noises are assumed to contaminate inlet 
duct pressure measurements. The cost function chosen to be minimized is a weighted 
sum of averaged square values of the output (throat exit static pressure) and its deriva- 
tive plus the control (inlet bypass door area) and its derivative. Optimal parameters 
were calculated for both proportional -plus -integral (PI) and proportional -plus -integral- 
plus -derivative (PID) single -loop controllers for a NASA designed mixed compression 
inlet. Feedback signals used were either throat exit static or diffuser exit static p es- 
sures. Designs were evaluated on the basis of averaged square values of output and 
control signals. PI control on throat exit static pressure proved to be the most effective 
compromise between complexity and ability to attenuate disturbances. Appendixes ai o 
included to describe the computer program as well as outline the solution to a sample 
design problem. 


A supersonic aircraft propulsion system consists of a supersonic inlet and either a 
turbofan or turbojet engine. The function of the inlet is to convert high velocity, low 
pressure air ahead of the inlet to low velocity, high pressure air more suitable for the 
engine's compressor. An efficient inlet for flight at Mach numbers above about 2 is the 
mixed compression type, having a convergent supersonic region followed by a divergent 
subsonic region. An increase in static pressure occurs in both regions as the flow is 
decelerated. A terminal shock separates the supersonic and subsonic flows. To maxi- 
mize efficiency, the terminal shock should be located near the throat in the divergent 
duct. If the terminal shock moves upstream into the convergent region, it jumps for- 
ward to form a strong shock wave ahead of the inlet. This occurrence is known as an 
inlet unstart, and results in increased drag and a rapid loss in pressure recovery. This 
in turn may lead to a compressor stall and/or combustor flameout. Conversely, if the 
shock moves too far downstream, the pressure recovery is reduced and distortion is in- 
creased at the compressor lace, which may also cause compressor stall. Thus, con- 
trols are needed to maintain the terminal shock close to but yet downstream of the throat 
for good pressure recovery without inlet unstarts. 

Figure 1 is a schematic diagram of a mixed compression inlet with a typical termi- 
nal shock control loop shown. The terminal shock is positioned using me bypass door 
loop. The door opens and closes to maintain a match between the inlet and engine air- 
flow as engine airflow demand changes. This tends to prevent the terminal shock from 
moving too far rearward, causing increased distortion, or too far forward, causing an 
unstart. Throat exit and diffuser exit static pressure signals shown give indications of 
shock position. They are fed back through the controller to drive the bypass door servo. 

The controller transfer functions are significant factors in the dynamic regulation of 
shock position. The design of such transfer functions has been studied in references 1 
to 4. A root locus design technique was reported in reference 1. A stochastic optimal 
control theory approach was reported in references 2 and 3, where the expected fre- 
quency of unstarts was minimized. A parameter optimization approach, where the pa- 
rameters in a fixed -form controller were selected so that the response approximated a 
desired closed-loop transfer function, was reported in reference 4. 

The approach taken in this report is also based on optimizing the parameters in a 
fixed-form controller. However, the cost function chosen to be minimized is the 
weighted sum of average integral square errors, due to a deterministic disturbance, and 
the mean square errors due to random noise on the measurements. This approach was 
motivated by the inlet control problem (ref. 5), where the compressor face disturbance 
is most conveniently described as a deterministic signal but where the noise on duct 
pressure measurements is definitely random in nature. The variables included in the 
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spective first derivatives. 

While the problem is formulated in the time domain, the optimization is done in the 
frequency domain. This is because system, disturbance, and noise models are most 
often obtained in frequency domain form. The conjugate gradient search (ref. 6) is used 
for the optimization. However, the standard procedure is modified to reduce search 
convergence time and insure that the resulting controller is stable. This is done by 
scaling and constraining the search parameters. 

In the next section, the general control problem investigated is defined and the cost 
function and gradient calculations are described. Then the parameter optimization pro- 
gram used is described, followed by a description of the inlet model to which the con- 
troller design method was applied. Finally, the results are presented, followed by ap- 
pendixes, which include one on the use of the computer program and one in which a sam- 
ple problem is presented. 


GENERAL CONTROL PROBLEM DESCRIPTION 
Description of Plant, Noise, and Disturbance 

Motivated by the supersonic inlet control problem, a control system structure was 
selected as a framework for the controller design problem. This structure is shown in 
figure 2. The linear plant (blocks G l and G 2 ) is assumed to have two measurements 
and x 2 available. (Symbols are defined in appendix A. ) Variable Xj, the outer - 
loop measureme.it. Is considered the output. The control is to consist of two blocks, 

Hj and H 2 , whose forms are specified at the outset, but whose parameters are to be 
optimized. The two inputs to the controller are assumed contaminated with Independent 
Gaussian noises Vj and v 2 , each having the same power spectral density $ v (u.'). The 
controller output drives actuator G^ to produce control u. The plant is acted upon by 
deterministic disturbance d which occurs at time zero. 


Cost Function Evaluation 


It is assumed that the purpose of the controller is to keep x^ and x^ as close to 
zero as possible but at the same time limiting the excursions in control u and control 
rate u. Thus, a cost function must be defined which adequately reflects the average de- 
viations in output and control caused by random measurement noises Vj and v 2 and 
deterministic disturbance d. We now proceed to develop such a cost function. 
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For convenience, define components of Xj and u such that 


x l = x lv + x ld 

x l = x lv + x id 
u = u v + u d 

u^v + ^d 




> 


( 1 ) 


where the first quantities on the right sides of equation (1) are the components due to 
measurement noise and Vg and the second quantities on the right sides are compo- 
nents due to deterministic disturbance d. First, separate cost functions will be derived 
for the measurement noise and deterministic inputs. Then, these separate cost functions 
will be combined into a single cost function. Consider first the case where only mea- 
surement noise is present. Define the cost function 


= lim — I q.x 

2 tJ t L 1 


l v (t) + l2 x lv^ 


+ r i u v (0 + r 2 


4*0 


dt 


( 2 ) 


The q’s and r’s are arbitrarily selected scalar penalties. Thus, C y is simply a 
weighted sum of mean square values. Next, consider the case where only a determinis- 
tic disturbance d is present. Disturbance d is assumed to be zero for t < 0. Define 
the cost function 


= T f [ q l X ld (t ) + q 2 x ld (t) + r l u d (t ) + r 2 u d (t) ] dt 


( 3 ) 


Time T is defined as the (arbitrarily selected) period during which control is to be 

effective in minimizing the effects of d. Factor 1/T is included to make C_, an 

d 

average integral squared quantity, and thus comparable to mean square quantity C 
Now, the total cost function can be defined as 


C = C v + C d 


( 4 ) 


Cost function C is to be minimized by proper selection of controller transfer function 
parameters. The problem may be simplified if we assume that T is large enough so 
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that x^, x^, u d , and ii d all go to zero before t = T. Then the total cost can be 
written as 



where the upper limit on the second integral C d«= has now been made equal to °°. 

Our aim is now to express C in the frequency domain, since, as was indicated pre- 
viously, we desire to perform the parameter optimization using frequency domain sys- 
tem models and data. In the following development, an expression for C will be ob- 
tained which is an integral in the frequency domain, over a range of 0 to <». The inte- 
grand will be obtained in terms of the following known quantities: noise power spectral 
density (PSD), plant, and controller transfer functions, and the Fourier transform of the 
disturbance. 

As previously noted, cost function C y (eq. (2)) is a weighted sum of mean square 
values of stationary random variables x ly , x ly , u y , and ii y . It can be shown (ref. 7) 
that the mean square value of any stationary random variable y can be expressed in the 
frequency domain as 

y 2 = lim — f y 2 (t)dt = — /%( w)d« (7) 

r— ■ « 2r J -t 2n J-o? y 

2 

where y is defined as the mean square value of y and $ («) is defined as the PSD 
of y. It can also be shown (ref. 7) that the PSD of the output of a single -input - single - 
output linear system whose transfer function is G(s) is given by 

* y («) = 4> z (w)|G(j«)| 2 (8) 

where 4> (u?) is the PSD of the input to the system. If the system of concern i3 a differ- 
z 2 

entiator, that is. G(s) = s, or y(s) = s z(s). then ®y(a<) - u? 4> z (u>). We can use this 

fact plus equation (7) to express C y in terms of PSD' s, obtaining 
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(9) 


C 


v 






dw 


where 4> x (w) is the PSD of x lv and $ u (w) is the PSD of u y . 

Next we must express in terms of frequency dependent variables. Given a 

time function y(r), having a Fourier transform y(ju>) f Parseval's theorem (ref,. 7) 
states that 





( 10 ) 


Then of equation (6) can be put in the form of t..e left side of equation (10) by allow- 

ing the lower limit in equation (6) be This can be done because the integrand of 
C doo is zero for time t < 0. If this fact is used in addition to y(ju) = jw y(jco), a fre- 
quency domain expression can be obtained for as 



Since the integrands of equations (9) and (11) are even functions of frequency, the lower 
limits can be set to zero, and an expression for cost function C obtained as 



To be able to numerically evaluate C, the PSD's and the absolute value squares 
must be expressed in terms of known system transfer functions, noise PSD 3> v (u.’)> and 
disturbance Fourier transform absolute value squared, |d(jo.')|^. Refer now to figure 2. 
First, variable x^ v , the component of Uj due to measurement noise, can be written as 
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x lv (s) = G a (s)G c/ (s)[h i (s)v 1 (s) + H 2 (s)v 2 (s)] 


(13) 


where 


G « (s)Gn(s) 

G d (s) F 

1 + G s (s)[g 1 (s)H 2 (s) + G 1 (s)G 2 (s)H 1 (sjJ 
Variable x^, the component of due to the disturbance, is 

x id (s) = G c; (s > d(s) 

Similarly, the component of u due to measurement noise is 


» = 


G c i 

G 1 (s)G 2 (s) 


jH 1 (s)v 1 (s) + H 2 (s)v 2 (s)] 


and the component of u due to the disturbance is 


u d (s) = G c ? (s)G a (s) 


’ H 2 (S) 

— + HJs) 


G 2 (s) 


d(s) 


Using equations (13) and (16), the fact that and v 2 are uncorrelated, and 
s=ja>. the PDS's of x ly and u y become 


<J> (w)= *>)|G a (jw)G c (jce)| 2 [lH 1 (ju)| 2 + |H 2 (Jw)| 

Xi_ r " L 


and 


<t> u (<e) = $ v (<e) 


G^ (j'i-’)G a (ice) 
G 1 (jw)G 2 (jw) 



Also, using equation (15) results in 

|x ld (jx ')! 2 - |G cZ Uo.)| 2 |d(jw )| 2 


(14) 

(15) 

(16) 

(17) 

letting 

(18) 

(19) 

(20) 
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and using equation (17) gives 


u d (M 2 = |G ci (]w)G a (jw)|‘ 


H 2 (ja') 


G^ja) 


+ H^jw) 




( 21 ) 


Substituting equations (18) to (21) into the cost function equation (12). cost C becomes 

( 22 ) 


c = C(b) = - J f F(-\b)du.’ 


where F(cc\b) is definet* as 


F(w,b) ft |G cZ U^b)i 2 ([Vq(cO + V r ( U )][|H 1 (i«.b)| 2 + |H 2 (i«,b)j 2 ] 


+ D («) + D r (u) 


H„(ju\b) 

-1— — + H.(jw,b) 


G 9 (jw) 


and 


V r (w) 


V q (u>) = (qj + a ’ 2 ^ 2 ) <I> v ( " ) l G a^ a) ^' 

A / 0 \ 1 G - (jw) 

= ( r i - «%> 


K,(") 


1 G 1 (jw)G 2 (iu)) 

1 

D q (w) = ^ (q L + w 2 q 2 j|d(jw )| 2 
D r (w) = - ( r i + u.' 2 r 2 )|d(ii’)| 2 |G a (ia>)| 




(23) 


(24) 


symbolism H.(j<,,b) and so forth has been introduced to indicate which terms are 
linn? nf controller parameter vector b. Note that V (tr). q a ’ r 


The 

functions of controller parameter vector b. * — -q' 

are all independent of the controller parameters. 
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Having now obtained an equation iron, which «o evaluate C as a function of known 
svstem pa-ametecs and t-a- fable controller parameters, the parameter optimiza.ion 
problem to be so, red is: minimize CfbJ 

r :r *:=c:idi — - 

■'ver the -roblem can be considered to be a completely stochastic one tf the determ.ms- 
;,c disturbance d(t) is rertaced by a Gaussian random variable with a o „ “ ; 

When » -at would be used in place of (1/T)!d(ial]- ,n equations (24). Also, he p, oo 

lem becomes entirely deterministic if measurement noise « is considered to 
In that case die noise PSD, * v (w) is set to zero in equations C 4). 


Cost Function Gradient Evaluation 

An expression for the cost funciioe. ieq. (22)) is not usually sufficient to allow ad 
effective solution to the parameter optimization problem. Most parameter op tmiza i 
, rrr-iHiPnK Fletcher-Povel!. steepest descents, etc. ) require also 
methods (conjuga e ' p 0 o( computing an approximate gradient is 

: rbrii: c » .he ^ ~ 

change in C divided by the chain in the respective b vector components However, 
a more efficient method is to derive ah explicit expression for the gradient XC m 
terms of b. Such an apprcach was alien in this report. Equation (22, can be used to 

write the gradient as 


VC(b) = - f FF(w,b)dw 
" n J 0 

Using equations (23) and (14). makes it possible to express VF(w,b) as 


(25) 



_ 

r 

x ■ 1 


VF(u\b) = 2 Real 

F(w.b)i 


^"•-’xH.Ow.b, 

G 2 (j«) ‘ . 




V. 


J 


^|G cZ (iu-',b)i 2 (D r (a.-) 


H 2 (ja',b) 


G 0 (jw) 


+ H l (jo.',b) 


H 2 (ja',b) 


G,(jw) 


h Hjtjw.b) 


+ |v cl (w) + V r (w)| • |H 2 (ja’,b)v|H2(iw.b)j 


+ HjUw.bjV 



(26) 
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where the fact that v|y(ju))j 2 = 2 Real [y(:a’)Vy(-ju;)] was used. Once the forms oi 
H (jit\b) and H 2 (jtL\b) have been chosen. :he required gradients can also be calculated. 
In\he computer program described in appendix B, Hj(s,b).and H 2 (s,b) are assumed *.0 
be of the following form: 


/ 

m l 2 / s 

H(3,b) = Kbs II — 

J =1 \ b 2i 


3/2 2sb^. 




Here, m m 2 , and rru are given. Exponents p 2 j an d P 3 j indicate whether me fa<_ - 
tors appear in the numerator or denominazor. are ± 1 . and are given for all j. ine 
parameter K could represent a transducer gain (for instance, conversion from a pres- 
sure measurement to a controller input voltage) and is given. With these restrictions, a 
closed form calculation can be made for VHj(jugb). For example, if H^s.b) has four 
parameters and is given as 



where = (b ^ , b 2 p b^^, then, 
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P rosjram m mg Considerations 


A computer program was written (see appeiutix B) which r u lutes the optimum l^t 
rumetors b using the conjugate gradient method (ref. G). In ~»ie program, the cost 
function (cq. (22)) and cost func i; on gradient (eq. (25)) are computed us ini? numerical in- 
tegration. In jxirticular. using the trapezoidal rule, the cost function can be u lit ton as 



2 - 





(30) 


where is the number of frequency points over which the integral is to be computed 

(S3), u.-. is the frequency at the i 1 ' 1 data point, +1 = u.'^ . and Wq - u'j. A similar 
1 d d 

expression can be written for the gradient. Judgment is required in selecting the spacing 

of the frequency points, such that integration errors are minimized. In addition, a var- 
iable step size feature is introduced to further reduce errors. This is done by using the 
fact that errors caused by a fixed step size are greatest near a system resonance. Reso- 
nance is defined herein as the portion where the closed-loop transfer function magnitude 
|G (s)] exceeds the open-loop transfer function magnitude. The program detects this 
condition and when it occurs. |G ( ,^ (ju.')|^ is evaluated at additional points within each 
frequency interval using interpolation between adjacent prespecified frequency points. 

This feature increases accuracy and also tends to prevent occurrence of another compu- 
tational problem, namely, that of the closed-loop system transfer function becoming un- 
stable during the search. Since the magnitude squared of an unstable G c ^ (j«) is the 
same as a stable one having the same pole magnitudes, the cost function calculation 
won't differentiate between the desired stable and an unwanted unstable solution. How- 
ever in becoming unstable during the search procedure. G c/ (jtr) will have an increasing- 
ly large resonant peak, which the program will tend to detect. Then the addition of the 
extra function values will ensure adequate accuracy in computing C and VC. That is, 
the decrease in system stability will be property reflected as an increase in the cost 
function. 

The search procedure used is similar to the conjugate gradient method of refer- 
ence 6. Two modifications were made to the conjugate gradient search to improve the 
speed of convergence in this application. The first was to constrain the search vector h 
components to not change sign. With the transfer functions defined by equation (27). this 
just means that only stable controllers are allowed as candidates for the optimal one. 

The second modification concerns scaling of the b vector. It is known that "spherical" 
cost functions (where the cost is more or less equally sensitive to each element of the 
b vector) tend to lend' themselves to rapidly converging searches. Thus, a scaling was 
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incorporated into llio search procedure sue!) that each parameter has approximately 
equal influence. This was done by defining a new scaled parameter vector p using the 
diagonal matrix A as follows: 


b = ,\ _1 b (31) 


where Aj is die magnitude of the i t!l element of b at die end of the previous intcration 
of the conjugate gradient algorithm. The cost function and gradient in terms of p are 
thus: 

C(b) = C(Ap) (32) 

and 

TF(p) = A VF(b) (33) 

Theory states that the unmodified conjugate gradient search will always converge for 
quadratic cost functions. In the modified conjugate gradient search, the coordinate sys- 
tem is changed at the start of each iteration. Thus in theory, the modified conjugate 
gradient algorithm may not always converge for quadratic cost functions. However, ex- 
perience to date on nonquadratic cost functions confirms that the aforementioned scaling 
algorithm gave convergence times less than or equal to those for nonscaled cases. Ap- 
pendix C shows the results for one test c;-se where scaling was particularly useful. 


15 = 


1/A . 


I.A., 


APPLICATION OF PAR AML :£R OPTIMIZATION TO INLET 
CONTROL DESIGN 

The plant, '.o which the parameter optimization method is applied in this study, is a 
NASA designed two-dimensional, mixed -compression inlet. A description of the inlet is 
given in references 8 and 5. Reference 8 gives experimental open-loop frequency re- 
sponses of the inlet's terminal shock and subsonic duct static pressures to overboard by- 
pass door area. Experimental frequency responses of the inlet with control are given in 
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reference 5. In these tests, one set of overboard bypass doors was used to generate the 
disturbance and a second set was used for control. 

Inlet and Noise Dynamics 

The structure studied in this report is shown in figure 3. The output Xj to be con- 
trolled is P 57 . a threat exit static pressure 57 centimeters from the cowl lip and down- 
stream of the terminal shock, which is used as an indication of shock position. Meas- 
ured signal x 2 is pressure P 8? . which is closer to the compressor face station, where 
the airflow rate disturbance d originates. Controllers Hj and H, drive bypass doors 
G which j\iss more or less flow to counteract the effects of dislutbunce d. Conti ol u 
is' u b>> . bypass door area (or airflow rate). 

The inlet dynamics and bypass door frequency responses were found in reference 8 
and are tabulated in table I. Magnitude data are shown normalized to the values at 
1 hertz. Experimental frequency response data were available over the range of 1 to 
150 hertz; points at 0.001. 300. and 600 hertz were extrapolated using transfer function 
models similar to those given in reference 8. 

From a limited amount of experimental data, measurement noises Vj and v 2 wore 
found to be uncer related and to have the same PSD >i> v (u>). The PSD is tabulated in 
table I, with the data for 300 and 600 hertz extrapolated to be equal to the value at 
150 hertz. 


Cost Function 

To simplify the discussion of results, the cost function of equation (12), by substi- 
tuting and collecting terms, can be written in terms of averaged square values as 

C = q^Xj + q 2 Xj + rjU 2 + u^IJ (34) 


where 


*i ■ ■ jf [%<“>♦?!*«<!“> 


(ov) + — |x. ,(jw)| 2 


dw = averaged square output (35) 


dw = averaged square output rate (36) 


13 


/*on 

»W 


u.') + 1 ; u d ( ju.'> j 2 


do.- = averaged square control 


(37) 


r 

" Jo 


*u u) + l - 


da' = averaged square control rate (38) 


Disturbance Model 


One of the critical assumptions in this study is die selection of the form of the dis- 
turbance d(s), The disturbance is assumed to have Uie form 


d(t) = Ae~ at ( 39 ) 

The disturbance repi esents a corrected flow rate change at the diffuser exit. For eon- 
vcnience, d(t) is taken as an equivalent diffuser exit area change (cm ) instead of cor- 
rected flow. A range of parameter a values is considered, from a = 4 to 400 radians 
per second, so as to account for our uncertainty in establishing the exact nature of the 
disturbance. 


RESULTS 


The majority of the results were obtained for a case designated as the reference 
case, defined as follows: 

(1) A single -loop control is used with H 2 = 0; H 1 is assumed to be of the following 
form: 


H^s) = K 




which is a proportional -plus -integral (PI) controller. PI control was chosen in view of 
the results obtained using PI control in previous inlet control programs. 

(2) An averaging time T of 1. 0 second is assumed. 

(3) The disturbance pole a is assumed to be 40 radians per second. Disturbance 
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ampins A is assumed to equal 8-1 square centimeters. This is about the same umplt 
tude that was used in unstart tests reported in reference 5. 


A v'c raided Square Value Com parisons for Reference Case 


Re-erence case designs are compared on the basis of their averaged square values. 
Three different cases of q 0 . r y and r, cost function penalties are considered. Since 
theva’ -sof q q 0 . r and r„ cannot all be varied independently, q, is set to one 
for all uses. For each ease, two of the penalties are set equal to zero and the third is 
varied from zero to infinity. Thus, although not all combination of penalties are c-am- 
ined, :re ones used will give a representative sample of all possible results. 

In figure 4(a). the normalized averaged square output rate (xj) N is plotted as a 

functio- of normalized averaged square output (x^) N - Quantities x“ and x“ are 
normalized with respect to their open-loop values. For each curve, the penalty indi- 
cated sees from zero to infinity in the direction of the arrow. The r { and r 2 cunes 
K o to the open -loop condition for rj or r 2 equal «. The q 2 = * case, which is 

equiva’ent to C = x?, appears to also have an open-loop solution, but in fact approaches 
it quite closely but doesn't reach it. Each of the three curves has a -’knee” type char- 
acteristic. thus a design trade-off exists between (*f) N and (i^) N .^For the q 2 case. 

for example, one might choose the design having (x 2 J N of 1. 2 and (x 2 ) N of 0. 5 as one 
which bis a fairly low value of (x 2 ) N while not having an excessively large value of 

(4_ T he q 2 curve lies below the r, and r 2 curves since it is the only case whore 
x 2 is penalized directly. 

In selecting a controller design suitable for actual implementation, it is important 
that the design not require control actuator capabilities beyond those available. Thus, 
control signal and control rate requirements are examined in figures 4(b) and (c)_Here. 

control and control rate are normalized to their values at the point of minimum x 2 (with 
r . r , = q, = 0). Figure 4(b) is a plot of normalized averaged square control as a func- 
tion of (^\ T . For constant (x 2 ) N . as would be expected, the case where u 2 is 


than the other two 


penalized Tj results in a controller which requires less ^u ^ 

cases. Similarly in figure 4(c). a plot ol against the r 2 curve tails be - 

low the r, or q 2 curves. Once the physical limits of control and control rate arc 
known e.g. , for the inlet, bypass door area, and bypass actuator power output) figures 


15 


m and (c) can bo used in selecting candidate controllers which would not require these 
physical limits to be exceeded. 

' Figure -1 can be used in the following_.itanner_to assist in coming up w-t.. a controller 
design "that minimizes a combination of and J* while not cau^ingvariables o 2 and 

U 2 to exceed limits. First, check to see whether the limits on x“. u . and u" are 
violated for the case q, = ^ - r, = 0. If the limits are not violated, tnen_this control- 
ler is acceptable . beingMhe one the assumed structure that minimizes x“ without re- 

'ard to variables x~. u 2 . and « 2 If one or more of the variables do exceed limits, 
determine the q,. r, and r, ralues at which the limits are reached ((q,),. and 

(r ) ) using figures 4(a). (b). and (c). respectively. Finally go back to the ..ompu e 

anVconduct a trial -and-error design using penalty combinations in the ranges 

and conduct , , <r until a design is found which minimizes 

(q 2 \ * q 2 " '• (r l } f “ r l ’ 1 2'/ “ 2 

x 2 and does not exceed the limits. 

' 1 The results displayed in [igure 4 can be examined In a conventional manner U dis- 
playing Unproportional and integral gains of .he reference case controller fo^vanous 

\ X 1/N’ 


values of lx 


This is done in figure 5. It can be noted that for constant 
proportional* 'gain is larger for the Fj case than for the r 2 case; 

gain bj is higher for the r 2 case than for the Fj case. That is. proportional gain 

m ° S Thre^rUwtor designs are compared on a normalized magnitude frequency re- 
sponse basis infigure 6. Each has a different type of penalty (q 2 . r y or r 2 ) but all are 

for a value of (x 2 lv = 0. 4. The frequency responses displayed are for G^G^jw) 

(open-loop response of Xj to d) and G^jaO (closed-loop response of ^ to d) It 

can be seen that the ^ design behaves more like a proportional controller in tha i 

attenuates the disturbance similarly at low and high frequency, *lule the redesign acts 
like an i.de^al controller since it attenuates the disturbance more strongly at .ow w.- 

qu^cies (less than 10 Hz) than at high. The case where x 2 is penalized In a compro- 

between the other two designs. 


tiffed of Disturbance Pole Location on A veraged Square Quantities 

The effect of haying a dislurbance which has a pole value a lar s er man or less than 
that of me reference case is shown in figure 7. For simplicity, comparisons are made 
only f or the „ cases with averaging time T = 1 second. It din be seen mar store the 

initial value of d(t) is me same [or all three cases, me open-loop value of x‘ is largest 
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for the smallest ralueof a. Also the case with the largest value of 
the largest value ct x^. as would be expected. 


d(tl (a = 400) has 


Ff fe^cf Adding Derivative Control Action and of Uei ngjir nor -Loop 
Cent rolle r on Averaged . Square Quantities 

In an attempt -a improve upon the single -outer -loop PI control, two additional con- 
figurations were investigated. Tin- first was the addition of derivative action to the 
outer -loop PI controller. The transfer function for this controller has the form 

Kb l (S/b >l + l )( s, b 22 + D Thc sccond was a „ inner-loop PI control, whirl: uses sig- 

nal x 0 , the diffuser exit static pressure. One reason for consider ing use of x 2 instead 
of x 2 is that x„ is nearer than x 1 to the point at which the disturbance enters. Thus 
it mi* T ht be expected that such a c« troller could better respond to a diffuser exit dis- 
turbance. In figure 8 both of these controllers are compared with the outer -loop PI con- 
troller for q 2 designsof die reference case. The outer-loop PI and PID controllers 

exhibit very similar (xj) N against (xf) N characteristics. Figure 8(a) shows the PID 
controller is able to reduce (xj)„ over that of die PI controller only at large values of 

(*i)n- Thc inncr ' l00p PI contro1 is not as effcctive in reducin S ( x i)n as either 01 thc 

olhcr^ controllers, except at low values of (if) N . Although there Is less phase lag la the 
loop between x, and d than between Xj and d (sec table 11; .here is also . leas 
Consequently the signal -to-noise ratio at the input to the inner-loop controller (HjW) 
is less than tor the outer-loop controller. It is believed that this poorer signal-lo-notse 
ratio accounts for the ineffectiveness ot the inner-loop control. It can he noted in fig- 
ures 8(b) and (c) that (?) N and (u 2 ) N are essentially identical tor the outer-loop PI 
and PID controllers. This, the added complexity ot the PID controller hardly seems_ 

warranted. Also. the inner-loop PI controller has poorer performance In terms of x, 
and X?. and also in terms of (u 2 )^ or ( IJ )n- 


SUMMARY OF RESULTS 

■ This report has demonstrated the use ot parameter optimisation techniques in the 
design of controllers tor a supersonic inlet. The basic problem formulation allows the 
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disturbance to be described as a deterministic signal but includes measurements which 
are corrupted with random noise. The controller design problem was set up as a param- 
eter optimization problem in the time domain but was solved in the frequency domain. 

A modified conjugate gradient algorithm was used to compute the optimum controller 
parameters. Control effectiveness was evaluated in terms of average square values of 
output, output rate, control, and control rate, and also in terms of frequency respon- 
ses. In applying the method to the inlet, it was found that, of the controllers investi- 
gated, proportional -plus -integral (PI) control using throat exit static pressure feedback 
was most effective. A proportional -plus -integral plus -derivative (PID) controller 
showed only marginal improvement over the PI control. A PI controller using diffuser 
exit, static pressure was inferior due to signal -to-noise problems. Also investigated was 
the effect or disturbance dynamic characteristics on controller performance. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 2, 1974, 

501-Z4. 
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APPENDIX A 


SYMBOLS 


A 

a 

b 

b l 

b 2i 

b 3j 

b 4i 

C 


d 

F 

G, 


disturbance amplitude, cm 
disturbance pole. rad sec 


'd<- 


controller parameter vector. ni p * l 


-m 


l 


M cl 


H, 


H, 


jw 


controller pain rad sec) 
controller polos or zeroes, rad, sec 
controller damping ratio 
controller natural frequency, rad/ sec 
total cost function 

cost function due to deterministic disturbance 

cost function duo 10 deterministic disturbance with upper limit set equai to - 

cost function due to measurement noise 

intermediate variable in cost function 

intermediate variable in cost function 
2 

disturbance, cm 


2 


cost functio' integrand or summand 

actuator transfer function, inlet bypass door actuator, cm'/V 
plant closed-loop transfer function, inlet throat exit static pressure to diffuser 
exit area disturbance, (N/cm )/cm 

plant transfer function, Inlet diffuser exit pressure fo diffuser exit disturbance. 

1 2 2 
N/cm /cm 

pl ant transfer function, inlet throat exit static pressure to diffuser exn static 
pressure, N/cm^/N/cm 
general controller transfer function 
outer -loop controller transfer function, V/N/cm 
inner -loop controller transfer function, V/N/cm 

integer, in eq. (27) 

Fourier transform variable, rad/sec 



K 



p 2j 

P 3j 

“l 

*2 

r l 

r 2 


T 


measurement gain, pressure transducer gain for inlet. V 'N/cm 

integer, number of elements in b vector 

integer exponent of free s's in controller transfer functions 

integer, number of first-order controller factors 

integer, number of second -order cont: oiler factors 

integer, number of frequency points in numerical integration 

2 

throat exit static pressme. N/cm 
diffuser exit static pressure, N/cm" 
transformed parameter vector. m p >'• 1 
integer. ±1 
integer, ±1 

penalty on output in cost function C 
penalty on output rate in cost function C 
penalty on control in cost function C 
penalty on control rate in cost function C 
Laplace variable, sec 1 

deterministic disturbance averaging time, sec 



time, sec 

2 

actuator output, inlet bypass door area, cm 
2 

bypass door area, cm 
component of u due to disturbance 
component of u due to measurement noise 
intermediate variable in cost function 
intermediate variable in cost function 

outer -loop measurement noise, noise on threat exit static pressure measure- 
2 

ment, N/cm 

inner-loop measurement noise, noise on diffuser exit static pressure measure 
2 

ment, N/cm 

2 

- plant output, inlet throat exit static pressure, N/cm 

plant inner-loop variable, inlet diffuser exit static pressure, N/cm 
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> . component of x^ due to disturbance 

Xj" component of x^ dee to measurement noise 

y dummy variable 

z dummy triable 

A diagonal scaling transformation matrix, m p x m p 
A clement of A 

r random noise averaging time, sec 

2 2 

power spectral density ot d. tom ) /Hz 

2 2 ' 

<5 power spectral density of u, (cm ) /rad/sec 

'ti ‘ ,22 

power spectral density of and (N/cin ) /Hz 

V 2 2 
* power spectral density of x^ v , (N/cm ) /Hz 

. power spectral density of y 

power spectral density of z 

o> frequency, rad/sec 

a?. frequency at i lh data point, rad/sec 

Subscripts: 

( )^ T averaged square value of, normalized to averaged square value for 
case with r^ = ^ = ( l 2 = a = *0, anc * T = 1 
( value at which limit is reached 

Superscripts: 

T matrix transpose 

-1 matrix inverse 

averaged value of 
derivative with respect to time 
* complex conjugate 


reference 
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APPENDIX B 


COMPUTER PROGRAMS 

This appendix describes the FORTRAN IV computer programs which mechanize the 
controller parameter optimization. The package consists of a main program, two sub- 
routines, and a block data subprogram, written for an IBM 360/67 TSS time sharing 
computer. The subroutine CCFM conducts the conjugate gradient search. The sub- 
routine CALFG computes the cost function and gradient. 


Dimensions 

The programs are dimensioned for a maximum number of frequency points N d of 
25 and a maximum number of controller parameters m p of 10. The vectors dimen- 
sioned N d are AD, Al, A12, DD, DQ, DR, GA, G1G2, HZ, PI, P12, PD, V. VQ, and 
VR, The vectors dimensioned m p are B, G, GS, ID, and Z, Vector H has dimension 

2m and vector W has dimension N . + 1. 

P a 


Main Program 

The main program handles data and performs preliminary calculations. The com- 
puter variables in the main program are defined in the Main Program Variable List. A 
flow chart for the mam program is presented in figure 9. The following is a FORTRAN 
listing for the main program: 


Main Program Listing 


C MA I N PROGRAM FOP CONTROLLER PAPM'FTFR C D T I M I ZAT 1 0‘! 

C0M°LEX G1(25),0102(25),GA(25) 

DIMENSION A1(25),A12(25),AD(25),P1(2S), D 12(25),PP(25),DP(25) 
DIMENSION U Z ( 25 ) , iff 26),V(25), v Q(25),VR(25),DO.(25),DR(25) 
DIMENSION 5(10), 0(10), 10(10), M ( 2 0 ) 

G0H , 0 , l/CALC/F0O,FV0,FnR,FVR,Gl,A102,GA,'',V0,VR,00,0R, 10,- 

1 K1,K2,MH1,NH2, I HST, NPATA, KMT, RN 
COW'ON/FMO/KOU'1T,KO 

COMMON /BLOCK /HZ, A1 , A12, AO, PI , P12, PO,V, A, AA,T,OAl , 0A12, GAP, OK, MP 
LOG I PAL KG 
EXTERNAL CALFG 

N AM C L I r T /NAf'3/A,AA,Al,A12,An,HZ,OAl,GA]2,GAP,GK,KPR,- 
1N0,P1,P12,P0,T,V 

NAMELIST /NAM2/ K1 , K2, NH1, MH2, B, in,LlMIT,Pl,P2,Rl, r> 2 
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c 
9 g: 


c 

9 C 

93 




' 4 ; 


0 O|MT DRH^R AM POP. R fC FPE rF 

u • |NST=(i*SFARn^ 2 »NAl* 2 , 3 -NAI' 3,4 = * 0 l <T) ,/- 

s' . I FO»(n*Cf»‘IV, 1-MOT '•O'SV^-POPOP)' ) 

3 ftofi p T, READ A»n PR KIT NV 3 
< (R [ T r (f, 03 P ) 

FORI^TC NAf‘3? 1 ) 

RFAO(5,:;Ar3) 

ipf/'iD rn.l) l)RITF (8 . NAM 3 ) 
c “ONVEnri^-i FRW *"*« R'SP TO COLLEY NO, 

NOATft = ':P 
r.fi= 0 Ai 2 

• )lrE v -(rM-i X(n.,=>nn*. 0174533 )) 

0102(O’*0A12*A12(l)*PFT o (CP c ’LX(0., 3 12(n*.0 1 7U533)) 

VJ( ! 1 =" 2 ( I )*F. 2831854 

nl ss c « : i:ssKiz;ffi *. .«* . » - . » 

WCWATA*l)-W(vnATA> 
r opol'PT, REAR AND «n«T »I M, 2 
51 VJR I TF(6, 340 ) 

343 FORMAT ( ' M ,M 2 ? ' ) 

REAO( 5 , MAf '2 ) 

SJSlTii*) «!.«.«.«« ,ii*n r« 

, 2 , FURMATC . 51 './. 

r n^'SiRS."!; «»’** "*«.<« nc , 

no 2 io i«l,nn 

ViVfaVK ! ) * *2 

ga$q=p c m_(9a( i ) * rnN jhc^a ( i ))) 

VR( I )!v < I )*(RltMW*n2)4{sn/RF4L(ei'i2(l )*rONJ'l(«102( I ) )) 
nn( i )=""( i )«(oi*WJ* rt 2 ) 

21 n ORC 1 )-RP( I )*(Rl^'VJ*R 2 )«GA 3 a 

r,n TO 8 7 

C oRPM°T A‘l n REAR INST 
777 WRITE(6,250) 

250 FORMMC INST?') 

READ (5,42) INST 

42 FORHAT(ll) lllCT 

r PRINT FlWlf^VSsPPVCr MEAOINR 
• ‘ ul 55 roi”T(tx!”i'.iJx.'/»/'.«. , 0 Ee , .iK , "“ l " , ' 5 ' > 

21 S I 7 . E * . 1 

EPS- 1 . E -5 
KNT-0 




%,■ 


’❖> 

%k-: 


c 

230 
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87 
7 51 


Wfmr(P,151> (M.n,!-- I/O 
r . n r“ v*. 7 ( 1 7 = ' , .3 3 7 0 H 1 . h) 

^ 777 

r M n 


A 

AA 
AD 
A1 
A 12 
B 


CALFG 


DD 

DQ 

DR 

EPS 

F 

FDQ 

FDR 

FVQ 

FVR 

G 

GA 

GAD 

GA1 

GA12 

GK 


Main Program Variable List 

disturbance time constant a (input variable) 
disturbance pulse amplitude (input var:uble) 

G a i]a') normalized, vector (input variable) 
normalized, vector (input triable) 

* G 2 (jw)| normalized, vector (input triable) 

controller parameter b. vector (input variable). When inputting. those in 
outer -loop controller transfer function H ^ must precede those in inner - 
loop transfer function Hg. Damping ratios must precede their natural 
frequency terms. 

subroutine (declared external in main program) which computes the cost 
function and gradient 

:d(jcr)| 2 , vector 

D (w), vector 

D r (a'), vector 

parameter change defining search convergence, for example, 1CT 5 

F = C = FDQ + FDR + FVQ f FVR, cost function 

costs of averaged square output and output rate due to disturbance 

costs of averaged square control and control rate due to disturbance 

costs of averaged square output and output rate due to measurement r ,isc 

costs of averaged square control and control rate due to measurement noise 

cost Junction (scaled) gradient, vector 

G a (jcc’), vector 

gain normalizing AD (input variable) 
gain normalizing A1 (input variable) 
gain normalizing A12 (input variable) 
transducer gain K (input variable) 
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GN GN - G A 12 


G 1 

G1G2 

li 

HZ 

I 

ID 


IER 

INST 


ITER 

J 

KG 

KNT 

fcoUNT 

KPR 

K1 

K2 

LIMIT 

N 

ND 


vector 

G t (j^ ) ■ C 9 (ja), vector 
storage, vector 

frequency in hertz, vector (input variable? 
index cf element in vector, integer 

integer vector which identifies corresponding jurameior in B us 'u *vp. 
(input variable). 1 - zero. 2 = gain. 3 - pole. 4 - complex zero damm.m 
0 = complex zero natural frequency Q - complex pole dumping, t • 0 ::: 

. ple*N pole natural frequency. 

search convergence parameter. 0 = convergence :o EPS in LIMIT: l = v< 
vorgence to EPS in LIMIT not obtained; 2 = probable error occurred. 

branching instruction parameter (input variable). 3 - search lor optimum 
2 - return to NAM 2 ir melist; 3 - return to NAM3 namelist; 4 - prim 1'r 
quency , system closed-loop frequency response normalized by GN). an 
system Nyquist plot. 

set equal to LIMIT except for INST = 4 case when FTER = 0. 
index of element in vector 

logical rariable TRUE, means compute gradient 
count of cost function evaluations 
count of line search iterations 

if equal to 1 causes NAM3 variables to be printed input variable) 
exponent of free s in hi (input variable) 
exponent of free s in H2 (input variable) 
maximum number of iterations (input variable) 

NH1 + NH2 

numbei of frequency points over which integration is performed N^ 


NDATA NDATA - ND 

number of parameters of B in HI (input variable) 
XU2 number of parameters of B in H2 (input variable) 

.G (jcv) in degrees, vector (inpu ’c) - 


PD 


PI 

in degrees, vector (input variable) 

P 12 

..G , (joj) ■ G 9 (jw) in degrees, vector (input variable) 

Qi 

q^ (input variable) 

N 

O’ 

q 9 (input variable) 

R1 

(input variable) 

R2 

r 9 (input variable) 

SIZE 

parameter step size; for example, set to 0. 1 at start < 

T 

T (input variable) 

V 

vector (input variable) 

w 

Vq(uj), vector 

TR 

V r (o>), vector 

W 

radian frequency u\ vector 

ww 

frequency squared 


Program Input and Output 

The program starts by printing a heading referencing the namelist variables and 
variable codes. Then the program prompts for NAM3 namelist data. The namelist 
variables are entered according to the FORTRAN rules for namelist data. The NAM3 
variables are A, AA, Al, A12, AD, HZ, GA1, GA12, GAD, GK, KPR, ND, PI, P12, 
PD. T, and V. Since eight of the variables are vectors, the input could be lengthy. 
Thus, an alternative to entering NAM3 data at run time is to use the block data subpro- 
gram. The following is a listing of the subprogram for the inlet investigated in the 
report: 


BLOCK OATA 

COMMON/ BLOCK/ HZ, Al, A 12, AO, PI, PI 2, P0,V,A, AA,T,CAl,CA12,C,An,GK,N0 
REAL HZ (25)/. 00 1,1. ,3. ,7. ,10., IS. ,20. ,25.,^*'., A3. ,50., GO. ,7(1.,- 
1- SO., 00., 10 a., HO., 120., 130., 140., 150., 300., 3*fi00./ 

REAL Al(25)/1. .997, .927, .818, .732, .643, .565, .411, .359, .326,- 
l“‘ i31S, .387. ,465, .486, .46, .407, .326, .303, .209, .042,3*. 0003/ 

REAL A12(25)/l.,l., . 95, .774, .6 07, .671, . 62 9,. 602, . 53, . 437, .347, .31 3,- 
1 . 334, .352, .3 18,. 3 22, . 299, . 28 8 , . 27 , . 26 , . 237, . 012, 3* . 0013/ 

PEAL AD(25)/1. , 1 .,. 099,1. ,.996, .993, .997, .9 96,1. ,1.012, 1.016,1. 012,- 
1” 1.01,1.011,1 .005,1.024, . 05 4, .90 8, . 782, .68 0, .609, .051, 3*. 0063/ 

. REAL PI (2 5 >/0., -1., -8., -IS., -26 ., -39 ., -46 ., -51 ., -58 ., -63 ., -f 2., -55.,- 
1' -49.j-55.,-65.,- 7 4.,- n 4.,-llo.,-125.,-115.,-98.,-174.,S*-178./ 

: REAL P12(25)/0.,-2.,-13.,-24.,-31.,-35.,-44.,-57.,-62.,-75.,-87.,-90.,- 
1-92. , -105. , -11°. , -128. , -145. ,-156. ,-16 6 . ,-172. ,-1 /7., -251. ,3»-26 6./ 
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P r M. e'M2S)/n.,-2.,-5.,-in.,-n.,-18.,-2V,-27.,-3?.,-4?.,-51.,-50.- 
REM. ''(25) / .74F-5, u*2 .4E-5,U . c -5,2*5 . 7 r -5, 7. r-5,*.F-5,4 .4E-5,- 

7 7.r-5,3.<4F-5,4*4.F-5,3*3.3F-5,U.4 c -S,7.r.5,3*u.e-5/ 

HMA 7,V,AA/1.,4'V,«4./ 

hua , j f ' # 0Ai^.Ai2 # ^n,^K/23 / .nir, .«:s n 2 ,^^. •, .^b:/ 

END ' ‘ 


The NAM3 variables are printed if KPR is se* to one. Then the program prompts for 
NAM2 namelist data. The NAM2 variables are B, ID. Ki, K2. LIMIT, NH1, NH2, Ql, 
Q2, Rl t and R2. Then the NAM2 variables arc printed and a prompt for the INST varia- 
ble is printed. 

The INST variable is entered in II format. The INST code values are 1, 2, 3, 
and 4. Making INST = 1 causes a search tor the optimum parameters. After the search 
results are printed, another INST prompt is issued. Making INST = 2 returns the pro- 
gram to request NAM2. Making INST = 3 returns the program to request NAM3. Mak- 
ing INST = 4 causes the system frequency response to be printed. The frequency Is 
printed under HZ, The system closed -loop frequency response normalized magnitude 
is printed under /Y/ and phase in degrees under DEG. The real part of the system 
open-loop transfer function is printed under REAL and the imaginary part under IMAG. 
Then the zero iteration search results are printed, and another INST prompt is issued. 


Subroutine CALFGfN, B. F, G) 

The purpose of this subroutine is to compute the cost function and Its gradient. The 
CALFG program variables are defined in the CALFG Program Variable List. Thuse 
variables, in the common blocks and subroutine call are labeled the same as those in the 
main program and are not repeated again. 

The Data statement IN and RSNAT values are variables in the variable step size in- 
tegration. Smaller step sizes are taken near a system resonance; that is, when the ab- 
solute value squared of the denominator of the system closed-loop transfer function 
1 1. + OPEN] is less than RSNAT. The variable IN is the number of subintervals into 
which an interval is divided. Upon detecting a resonance, the program fills in extra 
points starting from the previous frequency by linear interpolation, using data from the 
present and past frequency points. A flow chart of the integration logic for the subrou- 
tine CALFG is given in figure 10. 
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FORTRAN Listing of Subroutine CA.LFG 


'*2021 I 


.70 

to 




P0» NT; Pf>CV I HE 
P! P r S0NANnE 


^t^.rot'T ! * ;c - r ur r (v^/ / r i) 
rc‘P! PX :(IO) # r .Al?5) / "l(2S),Clo:(n5) 
r nr^t rv Y r.c rT S, 0°PV'l ! , ^1^2 1 , t , 1 

ri I j 2 fO /7PC5) 25), n 0( 25) / n P(??) 

n|"f‘!^!^’J np ( ! 0 ) , 1 n ( 1 0 ) , A ( 1 ), r C 1 ) 

p pro*nM / r .\ ^ p / rno , F V 0 , rnn 7 , n 1 , 0 1 G? , r A , \ ! , VO , V P , n 0 , HP , , - 

1 *M , K7/JP1 /!M2 , l N.^T/'^ATA, <MT / r iN 
pp'7‘0 41 /n'f / KOI’ NT, X r » 
l_ op ] c a » ‘"',P*‘NT 

npr R5VAT, IN/.5, >/ 

*1171 'V. I ‘ T ! P*!S 
KNT = F %, 7 ♦ I 
PFNT^Fl’E. 
no 7 0' J*!/' 

0(J)-0... 

A., 0. ) 

M2 ■OMPt.X ( 0 . , H , ) 

rnn«o. . 

P'Q»fl . - 
Fno»n, ’ 

FVR*n. . 

. v*>wm 

pyM«i n 

n-l ,/ c t '.)AT( IN) 

; pinrr i fop F’%c Lf open too* oma 

c FOR I M T FR°OlAT I AM RFTW rp N °0INT 
Pp 200 1-1/fOATA 
90 |OX-l n XN 
100 Pi -FLOAT M OX) *0 
RIM-1. -Rl 

I -1>*R1M*W( I )'°\ 

>CM n LX(0., W! ) 

c OP I in rONTRCU r ) 0 AND W2 
I F ( N n l . OT , 0 ) 

* I F ( TIM 2 . OT. 0) H2«S**F2 
. . 00 10 J-1,N 
nj«in(j) 

00 Tn (1,3, 1,2, 10,2,10), nj 

1 sT-(n/o(J)+i. )**{2-i n J) 

I F ( KO ) Z(J)-0/CS + P(J) )-(Pl -^AT( 1 OJ5-2 . ) 
no TO 8 , , 

2 sT = (s/n(j*i))**2 + ?.*s*R(>n/fl(J + i-' + i- 

1 F ( . "DT . KC) r.O TO OR 

*S*R(J)/(B(J+l)*S T l*(3.-RLOM(inj)) 

98 ST»ST« * ( 5- I H J ) 

GO TO ? 

3 Sf = B(vl) 

7f.J1 =Pf^| V (J . , 0 . ) 

t 1 F( J . LF . “! M 1 ) H1-M1.ST 

I r ( J . . ‘J 1 * 1 ) M2»M2*ST 

"\ r criTl*!"F 

FOR or <* 0 " ■' "oc • ROrniiTr SYSTFt’ OOSTG 
oa ( I -1 )»R i f'+OA{ | 5 *R| 

nn ( i -l)«Rlfi + r l( : >*11 

2 I =G1G2( I - l)*RIM*flV32( I )*RI 


CMF r K 

r» I; 


Gl* 
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110 

120 


125 

130 
14 0 


OPEN-G* t*(ftl |*H2*0K.'»|*H1) 

IP(!0X LT IN) O n TO HO 

I p ( rca I ( ( 1 . ♦0P c N) *PO M JO( 1 . +OPEM ) ) . LT . RONAT ) 00 TO 120 
'■? ft l»W( !♦!) 


IOXN-IN 


RONT" . FA LGF , 

WPl-CVfO )*F| OAT ( in»l )*W( l-l)*FLO/'T(IN-l n X-l))*R 

00 TO 140 
|OYN*l 

I F ( R8 M T ) 00 TO 130 

1 F ( I . OT . 2 5 W5S»W( 1-2) 

CR0T*-1 . ♦( (W( I ) ♦FI OAK 1 N-l )*V/( I -1) ) *R-WS$) / (W( I )~WSS) 


Fno»Fnn + 0P0T*P0V.' 
EVO*FVO+ f 'P. r T*VOV , *UMAO 
FnR*FnR + OROT*DR'. l *MU 
EVR*FVR«‘PR r 'T*VR 1, *UI , A r 
00 125 J-l/l 

0(J)»0(J)^R0T*0n(J) 

R S NT = .TRUE . 

00 TO 00 

WP1*W( 1 )*(1.-R)*W( I *1)*R 
Y»01G2l/(l.+0PEN) 

YY*R C AL( Y*0Om.J0(v) ) * (WPl-WS ) 

M2021 I »(^ll*H2)/0102t+Hl 
OQ'.MP.I *nr»( I )*R»M*D0( l-l))*YY 
VQW=(RI *VQ( I )+RIM*VO( I -1) )*YY 
VRU.(RI *VR( I )fRIU*VR( I -1) )*YY 
OR'.Mnt *PP.( I )+RIM*OR( I -1) )*YY 

HMAO-PPAL(M1*COMJO(H1)+H2*CONJO(H2)) 

MM»PFAl. ( M202 I I *C0*Ufi(M202 II)) 


*\ 



P5naF n o*nnv; 

FDR*FDR + PP. , !*UU 
FVO*P'/P-fVOl'*of'A.O 
FVR»FVR+VRW*Mt‘AO 
IF(.NOT.KO) 00 TO 180 
C COMPUTE ORA^I^UT 

5jo ( ^rvj+i'MA^* ( VO.' .'♦V r> ' /) + r)R'/* , ll') '■ OOMJO( -Y*GAI )+0RU*M202 I I 


ifcj c q.d r,F- ( rr* ( vnw»Rw> *H1 ) *COM JO ( HI ) 
IFCJ.EP.MH1 + 1) OF* ( 5T •OOMJO (01 I /0102 I )♦( VOW+VR'-') 
ns( Jl-P.EAl ( p * 00 ’Mn(Z(J) ) ) 


160 0(j)*0 (,))♦'' v!) 

180 |P( IMOT.LT. 4) GO TO 100 

0 pr | mj SYOTFU F n F.Q! , EM0Y RESPONSE ^NO OPEN LOOP 


MZ-\M/6. 2831054 
YfAOaC ARG{ Y/ON) 

nFO*ATAI'2 (AIM A 0 ( Y) ,RFAL(Y) )*57. 20578 
HR I TP (8, 300) HZ , YMAG, nc G, 0°E‘>! 


ion f nx- i nx*i 


*H2)*OOMJO(M2) 


W5*'M 

|C( inx.LE. !•!) 00 T" 100 
200 CONTINUE 

F-PDO + FV0-*-FnR*F''P 
RETURN 

300 FORMAT (1PE10.3,2(5X,2E10.3)) 


ENO 


If (tic user desired lo use a search routine other than subroutine CGFM, the scaled 
gradient AVC(b) may not be desired. In such a case the unsealed gradient may be ob- 
tained by dividing eacli term in the gradient vector by its corresponding parameter in the 
b vector. 


DEG 

DQW 

DRW 

GF 

GS 

G1I 

G1G2I 

HM 

HMAG 

H2G2II 

HZ 

HI 

H2 

I 

IDJ 

EDX 

IDXN 

IN 

J 

OPEN 
R,RI, RIM 
RSNAT 


CALFG Program Variable List 

ZG c ^ (ju?) in degrees 
partial product 
partial product 
gradient of F 

saved gradient partial sum, vector 

G-^jw) value 

GjCjw) • G 2 (Jw) value 

temporary value 

temporary value 

temporary value H2/G2 + HI 

frequency value 

H 1 

h 2 

frequency index 
ID(J) value 

counter for data insertions 
counter value for logic 

number of subintervals inserted into an interval near a system resonance; 
for example, IN = 3 means each interval is divided into thirds 

parameter index 
G a ■ GjIHj + GjHj) 
temporary values 

n 

condition defining existence of system resonance for 1 1 +OPENj '""RSNAT, 
that is, for RSNAT = 0. 5. The closed-loop frequency response magnitude 
is >1. 42 = (l/y'O. 5) times its open -loop-value. 
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RSNT 


temporary logic; true - resonance 


s 

s 

ST 

temporary value 

VQW 

intermediate product 

VRW 

intermediate product 

\VI 

u.' value 

\VP 1 

next frequency 

\V5 

saved WI 

wss 

saved WS 

V 

g c , ( k> 

YMAG 

|G c/ (ju.-)!/GA12 

YY 

partial product 


Z partial product in controller gradient, vector 

Subroutine CGFM(CALFG. N. B. F.G. SIZE.EPS. ITER. IER. 11) 

The purpose of Subroutine CGFM is to perform the conjugate gradient search func- 
tion minimization. Several nonstandard modifications relative to the conjugate gradient 
search described in reference 6 exist in the CGFM subroutine. The gradient is the 
scaled gradient .WC(b) and every iteration updated b parameters change the scaled 
coordinate system. This simplifies calculations somewhat but violates thcoictical con- 
vergence arguments for quadratic sosts. However, away fi'om the minimum point, 

C(b) may be quite nonquadratic; and around the minimum, changes in scaling are gener- 
ally small. This, it is believed, accounts for the decreased convergence times obtained 
using the scaling modification. Another nonstandard modification to the search is that 
the signs of the parameters b are not allowed to change during the search. This pre- 
vents formation of an unstable controller during the search. 

The CGFM program variables are defined in the CGFM Program Variable List. 
Variables carried over in the common block and subroutine call are labeled the same as 
those described in the mqfn program and are net repeated again. Figure 11 is a flow 
chart of CGFM. 
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FORTRAN Listing of Subroutine CGFM 


r.iiPnnuT | Mf rnrf'( r A I. F ^ , M , 0, F,0, S I 7.E, EPS, ITER, t ER,H) 
n|Me*'S|' v ' »»(].), u( 

L0nir«L KC 

Onf'MO'i/P'O/KOUNT, KO 
C Iti Tl A| ! 7ATI n ‘ ,t ' 

K0nr.'T=n 

CJCDr 2 . 

ifp »-5 
5 OFTA-O . 

NCYC-0 

is K*n 

no 20 

20 M(J)-DC.I) 

KO-.TRDF. 

CALL CM Pft(N,R,F,rt) 

C TEST FO n ST noD l '0 S p APCH 

I F ( KOI’NT . OF . I TFR) I FR= 1 
I F ( S I ZFI . LT. EPS) I C R- 0 
40 I F ( I FR.ftT.-2) RETURN 
c Cnrpt'TF OAST np/n ! FNT '.'EIGHT INS 
TSOR-O . 

00 50 J = l,*l 

50 TROR»TSPR+ r (»0**2 

iFtN'Ovr.co.n) GO TO 60 
BETA-TSOR/TSAVE 
60 S C A !_ r - 0 . 

DO 70 J-1,‘1 
J'! = J+N 

u (J‘0*-n(J)+P« r TA* l KJN) 

?o SCAi.r = o r Al r + A0S(H( J* 1 ) ) 

IF(S rA L r .OT.0. ) 00 TO SO 

1 c P-0 

00 Tn 40 

SO 50,MF®SIZ r /SCALF 
TSAVF-TSOP 
•|0''0* M 0 W 0 + 1 
r p o a r 
l-l 

0 U n DAT c P.'S 

100 00 110 J-1,H 
JN- J+ M 

110 P.( J)*"( J) * A OT, (1 . + S0ALF*"' J' 1 ) ) 

FS-F 

KO- . FALS C . 

0 All ''AL« v ‘(H,B,* # r ) 

C i_or | r r»!A #| nrc L'NF SFAPCM STFO f.lZF OR CO'.'CLUOFS SEARCH 
IF(K.OT.O) 00 TO 120 
IF(F.I.T.FS) r O TO 130 
IF(L.OT.l) 00 TO 140 
00 TO 150 

120 IF(F.LT.FSS) 00 TO 140 
GO TO 160 
130 SCAT/ » SCALE 
L- L + l 

SCA! c = SCA | F *otf n 
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FSS-FS 

IFCL.LT. 15) no TO 100 
I ER = 2 
GO TO 40 

C FIT OtJAHRATir CURVE TO 3 PTS . BRACKETINO LINE SEARCH f’IN. 
140 DO 148 J ■ 1 , N 
JM-J+M 
R1»H( J) 

R2*H( J )*ADS(1 .♦SCASV*H(JM) > 

R3»B(J) 

IFU.OT.3) R1*R1*(R2-R1)/STEP 
X1»(FSS-FS)*(R1-R3) 

X2«(FSS-F)*(R2-R1) 

|F(ABS(R2-R3).OT. EPS/4.) GO TO 147 
IF(L.OT.l) B( J ) =R2 
GO TO 148 

147 B(J)-(X1*(R1+R3)+X2*CR1+R2))/((X1*X2)*2.) 

148 IF(B(J)*H(J).LE.O.) B( J ) ■- . 1*B( J ) +EPS*R3 
C UPnATE SEARCH VARIABLE? 

SI ZE-SI7.E*(FLOAT(L)+2.)/4, 

KOlJMT-KOUNT+1 
I F (NCYC.GT.N) GO TO 5 
GO TO 15 

150 SCASV" SCALE 

SCALE* SOALE/(l.+STEP) 

K° K+l 

S ! Z Fb SIZF/(1.+PTEP) 

no to ioo 

ISO S I 7 Fa SIZF/C1.+STEP) 

DO 180 J**1,M 
180 B(J)-M(J) 

GO TO 5 

END 


CGFM Program Variable List 


BETA 

FS 

FSS 

J 

JN 

K 

L 

NCYC 

R1,R2,R3 

SCALE 


conjugate direction weighting 
saved F 
saved FS 
parameter index 
J + N 

indicator for step size reductions 

number step size increases within iteration 

number of iterations before restarting conjugate search 

terms in quadratic curve fit 

step size scale factor 
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SCASV 


saved SCALE 


STEP 
TSAVE 
TSQR 
XI, X2 


step size 
saved TSQR 

squared gradient terms sum 
partial product 
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appendix c 




SAMPLE PROBLEM 


A sample problem, using the inlet data from table 1, Is presented. The following Is 
a computer terminal printout of the problem solution: 

,T,"1 


1 

2 

3 

4 

5 

6 

7* 

8 

9* 

10 * 

11 

12 

13 

14 
15* 
16 

17 

18 
19 


r:ST*(l»SFAR0'!, 2«W2,3-NAf*3 / 4«PR' M i ) 

i pr s (o*coiv,i“NOT coMv, 2 »FRS np ) 

»ia:*.3? 

&nan3 *iend 

!, Si2 kl--l,nhl- 2 ,b-?n # l«OflO # M- 2 | ,l 5 ltrj«t- 5 n / 

ql-l # q 2 -. 0001 ,rl- 0 # r 2 -n,k 2 -n / nh 2-0 AenH 

v i v o him »jm2 - L I M 1 7* *1 ^ '• n •* > 

el'.&M.M* ».«•* *0 ".I* 

b= i.oooof ni i.oooo c 04 

risT? 

F» 3.11901 510E-01 I CP - B 9 
F90,FDR,FVn,FVP= 2.P84F-01 0.000 

9 _ _ _ - - . m r r V An 


n.on 


KMT , KOI INI, SITE- 102 

2.357F-02 0.000 


| r> s 


74 i.3fne-06 


4.7216F 
I M .ST? 


01 1.3050F 02 


20* 

4 


21 

HZ 


22 

1.000F-03 

1 . 

23 

1.000F 00 

1. 

24 

3.000E 00 

5. 

25 

7.000F 00 

7 * 

26 

1.O00E 01 

7 . 

27 

1.500F 01 

7. 

28 

2.000F 01 

fi . 

29 

2.500E 01 

G . 

30 

3.000E 01 

5 • 

31 

4.000F 01 

4 

32 

5.000E 01 

3 

33 

fi.OOOP 01 

3 

34 

7 . 000F 01 

3 

35 

8.900E 01 

3 

36 

Q. 000 E 01 

3 

37 

1.O00E 02 


38 

1.100* 02 

7 

39 

1.700E 02 

7 

40 

I.300F 02 

0 

41 

1 . 4 00 r 02 

2 

42 

1.500F 02 

7 

43 

3.000E 02 

1 

44 

fi.OOOE 02 

1 

45 

F« 3. 1190151. OF- 

46 

Fnn,FnR,FVQ 

, F W P" 

47 

0- 4 . 72 IGF 

01 7. . 

48 

MT.T? 



/Y / 0F0 

F31F-04 R.090F 01 
P06F-01 7.3R4F 01 
111F-01 5 . L 7 7 C 01 
4P-9C-01 1.248F 01 
705P-01-R.2P.1F. 00 
023E-01-1.75OC 01 
, 73 2 r -01-7> . 94 7F 0], 

> R71F-01-4. r 'll p 01 
.PR4F-01-r..73 n f 01 
,°74F-01-7. ft 13F 01 
.S30F-01-P.4PO I: 01 
*. 4 25F-01-P.PPGF 01 
.080^-01-9. 1R5 C 01 
P7PP-01- 1 . n R7 c 02 
.41 or -01-1 .??2 C 02 

t 37Qs>0!-i.324F 02 

In75c.nl- 1,40 4 F 02 

* 7 n n p-01-1, F°r r 02 

‘ 5P r )P-01-1 .P70F 02 
. 4 P, 8 F - 0 1 - 1 . 7 2 0 F 02 
1 .083F-01-1 .771 c 02 
. 2 OOF -02 1 .090 c 02 


rfm I MA r 

2 . G 2 9 F. - 0 1 - 5 . 4 G 2 F 
-1 . 1 8 P R -fl 1 - 5 . 4 G 7 c 

-2.OR7F-01-1.720F 
-1.R00F-01-R.145F* 
-1 .''21F-01-3.005F- 
-P . S36F-02-? . P50 p - 
-9.300F-02-7.1P3F- 
-1 . 13 RP- 01 - 1 . 7 n 4 F- 
-1 , OROF-01-1. 322F- 
-1 . 066F-O1-7 . R1R C " 

-n.4F4F-02-^.339P 
-P .622 e -07. -l.P1 7C 
-9. 7 33 F -02-5. 5° 3 C - 
-0 ! 2G4 r - n 2 2 . 75 5 c - 
- 0 , »rf c -02 5.1R9F- 
_I,^p4c_ 02 7 . 3 1 R c - 
7.7ROC-03 7.R27F- 
o > <n pp -02 R . 3 ? 7 e - 
4*.373F-02 3.5'*1 c - 
4 .It PR c -0 2 1.504P- 
3 . 0 27 r -02 1.M2F- 
_o # 70 Rr-n 5-1 . 3P,3 r - 
- 1 *007P-OR-° . n 75 c - 
0 ' 


H | FR* 1 KHT , K0MIT # 9 1 7.F" 1 

2.8S4P-01 0.000 2 . 352 p -02 0.00* 


03 

00 

00 

01 

01 

01 

01 

01 

01 

07 . 

02 
02 
03 
02 
02 
02 
02 
9 7. 

02 

02 

03 

94 

07 

.OOoe-Ol 


*User input. 
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Responses in capital letters ate program output (43 lines), and lower ease le el s u 
user input (5 starred lines). The line mm, hers to the left of each line were added for dis- 
cussion purposes. The listing shows the problem of optimizin' a two-parameter control- 
ler for a given cost function. The controller was a single-loop proportional plus integral 

. . _ . * • / IlT I. ( /U I l\/c — i\\ / -A fll- /bm ). 


controller with the form H 2 (s) = 0 and H^sJ/K -- bj 


h , (.s/I) 


+ 1)/ •*■> * (bj/s) H fh 1 /^ 2 1^ ’ 


M — 

The cost function was C = x* + 0. 0001 x\. It penalizes system output and output rate 
hut not control or control rate. 

Lines 1 to 5 were printed by the computer after the user called the pro-ram. These 
lines list the NAM3 and NAM 2 variables and parameters ID, INST, and IER. At line 6, 
the computer prompted for NAM 3 namelist. In line 7 the user signed in and out of the 
namelist without updating any values from the block data subprogram. Since K J • 
there is no print of the NAM3 variables. In line 8 the computer prompted for NAM 2. 
lines 9 and 10 the user supplied NAM 2 input. The initial parameter estimates were 
b = 10 and b,,, = 1000. The transducer gain K is part of the NAM 3 input and not in- 
pit here. In lines 11 to 13, the computer outputed the updated NAM 2 values, and in 
line 14 it prompted for an INST input. 

In line 15 the user entered a 1 to search for the optimum controller parameters. 

Lines 16 to 18 list the search results. The cost function F was 0. 3119. IER = 0 
meant that the (scaled) parameters were changing less than EPS d.E-S) as tanbe vert- 
tied noting that SIZF. was 1. 39 E -b. The Iteration line search count KOUNT wm 24. and 
the number of cost function evaluations KNT was 102. The q, and q 2 costs due to the 
disturbance, FDQ, was 0.2884 and due to the noise, FVQ, was 0. 0235 FDR and FVR 
costs were zero because r t and r 2 were zero. The optimum parameters b, = 47 

and b 9 , = 131 were returned. 

At line 19 the computer prompted for another INST Input. The user input INST = 4 
in line 20, which directed the program to display the system frequency response. 

Line 21 contains the column headings for the data given in lines 22 to 44. The first fre- 
quency was made 0. 001 Instead of zero (in NAM3) because H,(s) would be infinite at 
zero. A system stability spot check can be performed using the Nyquist criteria and c 
open-loop REAL and WAG data by reading the WAG column until a sign change occurs 
Twh ich Uie REAL par. is negative (between WAG . -0. 0056 and *0. 028). A stability 
check is that the REAL part is greater than -1. 0: -0. 09 is. A complete test may no, al- 
ways be so simple but lor this problem the response sin.ply spirals into zero. 

Lines 45 to 47 repeal the list of costs and search parameters, and at lino 48 the 

computer prompted for another INST input. 

This example problem is also used to study convergence and parameter scaling. In 

figure 12 the cost function contours are plotted as functions of 1 >j and b 2r The con- 
tours form an elongated valley more sensitive to b* than .to b, and .thus not -1 
scaled away from the minimum. The search trajectory starting at (10. 10 000) * 1. 
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shown The search run In C Uluru 12 look 3. 5 seconds oonlrul processor unit Him- on the 
m M 300/07 TSS. The same problem was run .Kin,' a standard conjugate gradient search 
[or two different hut constant sealing rules. One rule sealed Hie parameters Icy their 
initial estimates (10, 10 000). The problem ran about lice same- time, hut it converged he 
a Slightly less accurate F . 0.3121 Instead of 0.3110. The other rule was to seam the 
parameters near the solution (40, 130). In this case tlieceardc could not advance sign. 
cantly from the initial point for b 12 which remained at 10 000 while b, moved to 40. 


} 

I 

! 
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TABLE 1. 


- MIXED COMPRESSION INLET AND BYPASS DOOR FREQUENCY 


RESPONSE DATA AND 














Normalized averaged square output. ( x 

dp Control js function oi output. «( * funtr ; f u?*ri of output, 

HyurH A. \nrr.ali/i‘<l averaged Mjiidn> output rji. r^ntr,.! and control ra1»A <jv functions of mr\ ^\M\ 


averaged square output , is q^, rj. a. «t r ? penalties varv from n t<< -»•; r eft* r*»n« *- uise 
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transfer Cost Gam, /pro, 
font lion fundion rod sec rad 1 sec 




Fifjure 7. - Comparison of normalized averted square o-.Ijm i and output r <* 1r h* 
controllers designed for different disturbance |ede h .'jtums; q? designs v»V 
Averacjmqtme T t second; H ? is» 0 ; Mj.s* prnjK.rlii-ivil plus rontrr! 
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No mailed eve raged square output rate. 



ia> Normalized averaged square output rale. <tp Normalized average square control. 



ic Normalized average square control rale. 


Figure h. comparison of normalized avr raged square quantities for inner and outer ic-up PI 
Controllers and outer loop fMD controller; q^ designs onls. Averaging ttnr • 1 second; 
dr»! ii fliiiP''*’ pole d 4 (j radians per second. 
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InitidM/.itions 

, 


Set KG true 
Call CAIFG 



Stopping criterion 

_/ 


Mo 



Compute gradient weighting j 


i L_, 

Update B*s 

— , 


Line search step size 



| Yes 


line search optimization 

1 

i 




Set KG false 
Cal! CALFG 


Figure 11. ■ Flow chart for subroutine CGFM. 
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